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Abstract 

We investigate the thermal activated magnetisation reversal in a single ferromagnetic nanoparticle with uniaxial 
anisotropy using Monte Carlo simulations. The aim of this work is to reproduce the reversal magnetisation by 
uniform rotation at very low temperature in the high energy barrier hypothesis, that is to realize the Neel-Brown 
model. For this purpose we have considered a simple cubic nanoparticle where each site is occupied by a classical 
Heisenberg spin. The Hamiltonian is the sum of an exchange interaction term, a single-ion anisotropy term and a 
Zeeman interaction term. Our numerical data of the thermal variation of the switching field are compared to an 
approximated expression and previous experimental results on Co nanoparticles. 
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1. Introduction 

The exciting field of nanostructurcd materials 
offers many challenging perspectives for fundamen- 
tal research and technological applications. Funda- 
mentally, the investigation of magnetic nanopar- 
ticles has already led to a better understanding 
of magnetic phenomena, such as domain struc- 
tures [1] and superparamagnetism [2]. From the 
technological point of view, their potential appli- 
cation in magnetic recording media [3,4] is respon- 
sible for a great interest. Indeed, sufficiently small 
ferromagnetic nanoparticles are single domain [5,6] 
and are good candidates for high density storage. 
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However, too small particles are superparamag- 
netic even at low temperature and no information 
can be stored. So an optimal size has to be found 
for each material. 

Moreover high density storage requires a precise 
knowledge of the magnetisation reversal dynamics. 
The investigation of the magnetisation switching of 
individual ferromagnetic nanoparticles driven by 
an applied field at very low temperature is now pos- 
sible using the micro-SQUID technique [7,8,9,10]. 
By waiting time, switching field and telegraph 
noise measurements on an ellipsoidal Co particle, 
these authors concluded that the magnetisation 
reversal can be described by thermal activation 
over a single-energy barrier as proposed by Neel 
and Brown [11,12]. On the other hand, because 
of their small size, nanoparticles are good candi- 
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dates for numerical investigations in which they 
can be studied with a realistic size. Magnetisation 
switching in small ferromagnetic nanoparticles 
has been previously investigated by Monte Carlo 
(MC) simulations on Ising [13,14,15] or Heisenberg 
models [16,17]. For highly anisotropic systems, it 
has been shown that the reversal process is not 
uniform but rather nucleation-likc. Analytical ar- 
guments and MC results have been compared with 
good agreement [14,15]. 

The simulation of the uniform rotation requires 
models with continuous degrees of freedom, like 
the classical Heisenberg model. Then, by varying 
the anisotropy strength or the particle size, it is 
possible to investigate the crossover from uniform 
rotation to nucleation [16]. In this paper, we wish 
to reproduce the reversal magnetisation by uni- 
form rotation in a single ferromagnetic nanoparti- 
cle with uniaxial anisotropy at very low tempera- 
ture. In our case where the spins are strongly cou- 
pled and the energy barrier is large compared to 
ksT, that is very large metastable lifetime, the 
Langevin dynamics formalism [12] would require 
a huge computational effort. Indeed, for such sys- 
tems, this method is restricted to timescales of the 
oredr of a few nanoseconds [18]. So we have cho- 
sen Monte Carlo simulations which do not take 
into account the precession process but allow to de- 
scribe considerably larger timescales. We have used 
the standard Metropolis algorithm. Unlike previ- 
ous works devoted to the lifetime of the metastable 
state [16,17], we are interested here in the switch- 
ing field which depends on the temperature but 
also on the field sweeping rate. Actually, we aim to 
check the scaling behaviour predicted by approxi- 
mated calculations [19,20] in the framework of the 
Neel-Brown model. Our results will be compared 
to recent experimental measurements [7,8] which 
are in good agreement with this scaling behaviour. 
We also wish to precise the limits of validity of this 
approximated expression in terms of temperature 
and field sweeping rate. The study is carried out for 
several particle sizes and anisotropy constant val- 
ues. In Section 2, we briefiy remind the main pre- 
dictions of the Neel-Brown theory and also derive 
results concerning the switching field. The model 
and the simulation technique are presented in Sec- 
tion 3. The results are discussed in Section 4 and 



a conclusion is given in Section 5. 



2. Magnetisation switching of a single 
nanoparticle 

The stable configurations of a ferromagnetic 
particle of radius R depend on the two ratios R/X 
and R/5 where A ^ \/j /M, is the exchange length 
and 6 ~ \/ J/D is the domain wall width (J, Mg 
and D are, respectively, the exchange interaction, 
the spontaneous magnetisation and the uniaxial 
anisotropy constant per atom). More precisely, 
the nanoparticle is single domain (with all spins 
aligned) if its radius is smaller than a critical 
value [5] . Such a nanoparticle exhibits two ground 
states of opposite magnetisation along the easy 
axis, for example the x axis. Under an applied 
field H in the opposite direction of the particle's 
magnetisation, the magnetisation is expected to 
reverse by uniform rotation, that is, the spins re- 
main parallel during the reversal. According to the 
Stoner-Wohlfarth model [21] and its generalisation 
[22] the energy barrier which is only due to the 
anisotropy of the system can be written, at K, as 



AE{H) = AEo 1 - 



H 
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where AEq, the anisotropy energy barrier in zero 
field of a AT-spin particle, is given by 



AEn = NDS"^ 



(2) 



(.S* is the spin value). The exponent a depends on 
the direction of the applied field: a is equal to 2 
when the field is along the anisotropy axis or per- 
pendicular to it and a « 3/2 else [7]. H^y^, the 
field for which the energy barrier vanishes, is the 
switching field at K . Assuming that the intensity 
of the atomic moment is related to the spin value 
by m = gusS leads to 
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If ksT « AE{H), only the up and down magne- 
tized states along the x axis can be observed and 
if fcsT << J, Eqs. (1), (2) and (3) remain roughly 
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valid. Prom the Neel-Brown theory [11,12], the Ufe- 
time of the metastable state can be expressed by a 
thermally activated expression 



t(T, H) = To cM^E{H)/kBT) 



(4) 



where the prefactor tq is usually supposed constant 
for simplicity [7]. Actually, it is the transfer of the 
thermal energy of the lattice to the system of cou- 
pled spins which is responsible for the rotation of 
the magnetic moments through the equatorial y— 
plane. So To depends on the temperature but also 
on the applied field [11,12,23,24]. 

The magnetisation reversal can be investigated 
by increasing the applied field over time and mea- 
suring the switching field. Since the probability 
P{t) that the magnetisation has not switched af- 
ter a time t decreases as the temperature increases 
{P{t) = ex.p[-t/T{T,H)]), the switching field de- 
creases as the temperature increases. On the other 
hand, the switching field increases with the field 
sweeping rate v = Since the magnetisation 
reversal by thermal activation is a stochastic pro- 
cess, the covered trajectory in the phase space is 
different from an experiment to another, so the 
switching field is varying and one has to deal with a 
switching field distribution. For small enough field 
sweeping rates, this distribution is given by [19] 



P{T,H) = 



t{T, H)v 



exp 



H 



dH' 



t{T, H')v' 



(5) 



It is the product of the probability ^(j^fj-^^ that the 
reversal occurs in the interval [H, H -f dH] by the 
probability that the reversal has not yet occured. 
Prom the Taylor's development o{p{H) around its 
maximum, one can deduce the mean switching field 

by 
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where c = ksHg^/ {aroAEo) and £ = 1 — 
Hsw/H^yj [19,20]. Then, assuming that the stan- 
dard deviation is roughly equal to the half-width 
of the distribution, one obtains [19,20] 
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This quantity increases with temperature and with 
the field sweeping rate. 



3. Model and simulation technique 
3.1. Model 

Our model consists of a simple cubic lattice in- 
cluded in a sphere of radius R with free bound- 
ary conditions. Each site is occupied by a classical 
Hciscnberg spin = (Sf , Sf , S^) of modulus S 
in order to reproduce uniform rotation and to in- 
vestigate the effect of the anisotropy strength. The 
hamiltonian H is given by 

H = -JY1 S,.S,-D^(5n'-H.^S, (8) 

<i,j> i i 

where J > is the ferromagnetic exchange inter- 
action limited to the nearest neighbors and D > 
is the anisotropy constant on each site. The first 
term favors alignment of all spins and the second 
one favors alignment along the x axis (up or down). 
The last sum represents the interaction with the 
applied field. Here, we have considered an applied 
field along the x axis: H = —Hex where e^, is the 
unit vector and H > 0. Then, the ground state 
corresponds to all spins being antiparallel to the 
unit vector e^ whereas the metastable state corre- 
sponds to all spins pointing up in the x direction. 
We have to mention that we have taken here the 
same exchange interactions, anisotropy constants 
and anisotropy axis for the core and the surface 
of the nanoparticle although surface anisotropy is 
usually considered as radial and stronger. Since our 
calculations are restricted to very small nanopar- 
ticles, much smaller than the critical size to be 
single-domain in the remanent state [25], we have 
neglected dipolar interactions. The magnetisation 
per spin of the particle is defined by 
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3.2. Simulation technique 

The numerical procedure is the importance- 
sampling MC method [26,27,28]. In our simu- 
lations, wc have used the standard Metropolis 
algorithm [29]. Each simulation is performed at 
a given temperature T. The particle is initially 
magnetized in the x direction by applying a strong 
field Ho = HoBx {Ho > 0) during n MC steps 
(a MC step (MCS) corresponds to the scan of 
all spins once, trying a rotation for each). Then, 
the antiparallel applied field is ramping by small 
jumps SH: Hp = —p SH Bx- Since the field is kept 
constant during n MCS, the field sweeping rate 
is given hy v = 6H/n. The simulation is stopped 
when the magnetisation reversal is observed (the 
arbitrary chosen criterion is < —0.85*) and the 
switching field value is stored. For each tempera- 
ture and each value of the field sweeping rate, a 
relatively large number of simulations Ug has to be 
performed to get a reliable estimate of the mean 
switching field Hsw{T,v). This has been done us- 
ing a parallel version of the code. The error for 
Hgyj[T,v) will be given by a{T,v)/ ^/nl + dH/ua 
where a{T, v) is the numerical estimate of the stan- 
dard deviation of the switching field distribution. 

Starting from the initial state magnetized in the 
X direction, a uniform rotation of the magnetisa- 
tion of an angle 6 produces, using Eq. (8), an en- 
ergy variation 

6E = NSil - cos6')[L>5'(l + C0S6') - H] 

where there is no exchange energy. Simulating this 
mechanism by a single spin rotation (SSR) algo- 
rithm implies an increase of the exchange energy 
at each individual rotation of 6. The corresponding 
energy variation of the system is 

SEssn = S{1 - cose) [zJS + DS{1 + cose) - H] 

where z is the coordination number. Since the ac- 
ceptance of the SSR requires SEssr < at JT, 
no individual rotation can occur at K if H < 
-^SSR ~ •Z'/S' -l- DS{1 + cosO) and, consequently, 
no uniform rotation of all spins. This clearly means 
that the algorithm fails a.t K since it gives a 
switching field value larger than the field which 
produces the flip of a surface spin {z = 1), that 
is H = JS which has nothing to do with the ex- 



pected value -ff°u)- very low temperature and 
small applied field as in our simulations (that is 
^-E'SSR > 0), SSR has a very small occurence prob- 
ability and uniform rotation itself can require a 
very large number of MCS. So, the switching field 
at very low temperature will be overestimated ex- 
cept for small enough field sweeping rates. There is 
a temperature range and field sweeping rate range 
for which the algorithm dynamics is suitable for 
the investigation of the uniform rotation. 



4. Numerical results 

The spin value has been fixed to 5' = 1 for all 
nanoparticles. The physical parameters N, J and 
D have to be chosen carefully in order to compare 
our results with previous experimental ones with 
a reasonable computational effort, that is small 
sizes. Actually, the magnetisation reversal process 
depends, as the stable states, on the two ratios 
R/X and R/S. In their experiments, W. Wernsdor- 
fer et al. [7] have investigated a single domain Co 
nanoparticle with a radius of (12.5±2.5) nm which 
yields R/X ~ 1.5 and R/5 ~ 0.7. Here, we have 
considered a 33— spin nanoparticle {R = 2a, a is 
the lattice parameter) with J — ks and D = 0.1 J 
so ii/A ~ 1 and R/5 ~ 0.3. These values are of 
the same magnitude as the experimental ones. Ac- 
tually, in our model, wc have miderestimated the 
exchange interaction of a factor 250 in compari- 
son to Co which decreases the two characteristic 
lengths A and 6 and allow us to investigate very 
small nanoparticles. The choice of J has only ef- 
fect on the simulated temperatures since we have 
imposed the low temperature condition ksT « 
J (0.01 < ksT/J < 0.1). It has to be mentioned 
that the quality factor Q — (X/S)^ which defines 
the hardness of the ferromagnetic material is equal 
to 0.09 in our model, of the same magnitude as 
the Co value 0.2. The mean switching field, at a 
given temperature and for a given field sweeping 
rate, was obtained by averaging over 100 different 
realisations. 

Firstly, we checked that our definition of the field 
sweeping rate v = SH/n is reliable; then, we have 
determined the thermal variation of the switching 
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Fig. 1. Thermal variation of the mean switching field for 
6H/n = 5 X 10"^ J and 5 X lO"'-* J for three different 
values of (n, 5H) (N = 33 spins and D = 0.1 J). The error 
bars are smaller than the symbols. 

Table 1 

Field sweeping rate versus the number of MCS for SH = 
5 X 10-3 J. 



n (xlO^) 


5 


10 


50 


100 


500 


1000 


5000 


10000 


100000 


v/J (xlO-'5) 


1 


0.5 


0.1 


0.05 


0.01 


0.005 


0.001 


0.0005 


0.00005 



field of a given particle with different values of SH 
and n such that 6H/n is kept constant. In Fig. 1, 
we have plotted the two cases SH/n = 5 x 10"^ J 
and SxlQ-'^J. It can be seen in both cases that the 
three curves corresponding to different numbers 
of MCS are almost identical indicating that 5H/ n 
is a reliable definition for the field sweeping rate. 
The fluctuations of the data are smaller for v = 
5 X 10~^ J. In the following, we have fixed 6H = 
5 X 10"^ J and only n is varying (Table 1). 

4.1. Temperature and field sweeping rate 
dependence of the reversal field 

We present here our numerical results con- 
cerning a nanoparticle oi N = 33 spin with an 
anisotropy constant D = 0.1 ■/. In Fig. 2, we have 
plotted the time evolution of the modulus and 
the components of the magnetisation per spin for 
ksT/J = 0.01 and 0.1. For fcsT/J = 0.01, the 
thermal fluctuations of the modulus are very small 
and the uniform rotation is observed (the modulus 
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Fig. 2. Time dependence in 10* MCS of the modulus 
and the components of the magnetisation per spin at 
ksT/J = 0.01 (left side) and fejjT/J = 0.1 (right side) 
{N = 33 and D = 0.1 J). For reason of clarity, the 

component is not shown on the right side. Note that we 
have plotted the time evolution for only three values of the 
field. 



is roughly constant during the reversal). The re- 
versal duration is about 4 x 10^ MCS which means 
that the reversal dynamics is slow. For UbT/ J = 
0.1, the reversal is still uniform but thermal fluc- 
tuations are visible. Many reversal attempts can 
be seen before the reversal itself whose duration is 
significantly reduced to about 4 x 10^ MCS. We 
can conclude that these parameters allow to be 
in the uniform rotation regime unlike the case of 
larger anisotropy constant, as it will be seen later. 

The probability distribution of the switch- 
ing field for 5000 samples is plotted in Fig. 3 
at ksT/J = 0.01 and 0.1 for two different field 
sweeping rates v/J= 10"^ and 10"^. We clearly 
see; the; gaussian shape of the distribution whose 
width proportional to a increases with v and T in 
agreement with Eq. (7) . 

The field sweeping rate dependence of the mean 
switching field for several temperatures is shown 
in Fig. 4. An almost logarithmic dependence at 
small sweeping rate can be observed as in Ref. 7 
whose experimental results are shown for compari- 
son. In these experiments the mean switching fleld 
has been measured over several hundred cycles. 

To check the validity of Eq. (6), one can 
plot the mean switching field values versus 
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Fig. 3. Probability distribution of the switching field at 
fesT/J = 0.01 and 0.1 for v/J = 10"^ and 10"'' (AT = 33 
and D = 0.1 J). Each histogram is obtained from 5000 
samples. 
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Fig. 4. Field sweeping rate dependence of the mean switch- 
ing field of a 33-spin particle with D = 0.1 J for different 
temperatures 0.01 < kBT/J < 0.1 (increasing temperature 
from the top to the bottom). Insert: Experimental results 
of Wernsdorfer et al [7] . 

1/2 

[fcBnn(^)] where e depends on Hgw If 
the scahng behaviour is satisfied, all data points 
should collapse on a single straight line by choos- 
ing correctly the constants and c. Then, a 
self-consistent fit is obtained when the vertical 
axis intercept is the chosen value Hg^^^. Prom Eq. 
(6) , the slope is -H^^/ (AEo)'^/'^ from which it is 
possible to extract an estimate of AEq. Putting 
back H^^ and AEq in the expression of c provides 



Fig. 5. Scaling plot of the mean switching field values of a 
33-spin particle with D = 0.1 J for 5x10-" < v/J < 10"^ 
and 0.01 < ksT/J < 0.1. The dashed line shows the mean 
linear fit. Insert: experimental results of Wernsdorfer et 
al [7]. 

a numerical value for tq. Such a fit is shown in 
Fig. 5. 

Only the curves corresponding to the smallest 
field sweeping rates {v/J < 10~^) collapse on a 
single straight line. Numerical data associated to 
higher sweeping rates move away from the scal- 
ing behaviour as the temperature decreases. This 
deviation is attributed to the inefficiency of the 
SSR procedure at very low temperature (see Sec- 
tion 3.2) because our numerical switching field val- 
ues can be larger than H^^ which is not physical. 
No deviation due to thermal fluctuations has been 
noticed even when ksT/J 0.1. The results of 
the best fit obtained with c/ks = 10~^ and H^^ = 
0.215 J are reported in TABLE 2. The vertical axis 
intercept is consistent with the value used to per- 
form the fit but is slightly higher than the expected 
one H'^.^ ~ 0.2 J (Eq. (3) with g^oiiB = !)• Our 
estimate of the energy barrier in zero field 2.8 J is 
lower than the expected value AEq — 3.3 J (Eq. 
(2)). 

4.2. Influence of the particle size on the reversal 
field 

Since the magnetisation reversal mechanism is 
very sensitive to the system size, we have studied 
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Table 2 

Parameters of the fits and estimates of tlic switciiiiig field 
at -ft' and the energy barrier in zero field for two particle 
sizes and 13 = 0.1 J. The axis intercept corresponds to 
H'^sw/J and the slope to -{H°^/J)/{AEo/jy/'^. 



N 


33 


123 




0.215 


0.215 


c/kB 


10-3 


10-4 


axis intercept 


0.219 


0.215 


slope 


-0.132 


-0.080 


AEo/J 


2.8 


7.2 


TO 


40 


150 



the influence of tlie number of spins on the mean 
switching field variation in the case D = 0.1 J. 

For N = 7 {R = a), the time dependence of 
the modulus and the components of the magneti- 
sation exhibits no significant difference from the 
case N = 33 except an increase of the thermal fluc- 
tuations. Indeed, for a given temperature, thermal 
fluctuations increase as the size of the system de- 
creases. A consequence is that the reversal duration 
is reduced to about 6 x lO'"' MCS and 10^ MCS at 
ksT/J = 0.01 and 0.1 respectively. For N = 123 
{R = 3a), the two ratios R/\ ~ 1.5 and R/S ~ 0.4 
are closer to the experimental values; in that case, 
we have observed again that the reversal is uniform 
which agrees with experimental results [7]. 

The plot of the mean switching field values ver- 
sus [ksTln (^)] is drawn in Fig. 6 for iV = 7, 
33 and 123. 

For N = 123, using the same procedure as in 
the previous section, we observed similar results as 
for A/' = 33. A self-consistent fit can be obtained 
with the same value of H^^ as for A^ = 33 in agree- 
ment with the fact that should not depend on 
the particle size (TABLE 2). The extracted energy 
barrier in zero field is significantly lower than the 
expected value AEq ~ 12.3 J again. 

In the case N = 7, the curvature observed as 
the temperature increases corresponds to the tran- 
sition towards the superparamagnetic regime. In- 
deed, the particle is still strongly magnetized at 
these temperatures (the time average of the mod- 
ulus of M at ksT/J = 0.1 is 0.9345), so the zero 
mean switching field (actually 6H) obtained for 
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Fig. 6. Scaling plot of the mean switching field val- 
ues of a 7-spin, a 33-spin and a 123-spin particle with 
D = 0.1 J and 0.01 < fcsT/J < 0.1. The values of 
the field sweeping rate are 5 x 10~® < v/J < 10~^ for 
AT = 7, 5 X 10-" < v/J < IQ-"^ for AT = 33 and 
5 X 10-1° < v/J < 10-"^ for AT = 123. 

the smallest field sweeping rates means that the 
magnetisation spontaneously reverses in zero field. 
At low temperature, the data do not lie on a sin- 
gle straight line. The deviation observed for v/J = 
5 X 10~^ and 10~® is attributed to tmlikely sin- 
gle flips of spins near the particle boundary (the 
coordinence number is only 1 for 6 spins) which 
induce the magnetisation reversal earlier as ex- 
pected. For this reason, we did not perform simu- 
lations with smaller field sweeping rates. This de- 
viation has not been observed for A' = 33 and 123 
because surface effects are much less pronounced. 

4.3. Influence of the anisotropy constant on the 
reversal field 

The ratio D/J is also supposed to play an im- 
portant role in the magnetisation reversal process. 
Here, we have investigated for a 33-spin nanoparti- 
cle the effect of increasing the anisotropy strength 
{D = 0.2 J, 0.3 J and J) to evidence the transi- 
tion to a non uniform rotation regime and to study 
its consequence on the switching field variation. 

For D = 0.3 J {R/6 ~ 0.5), the reversal is not 
perfectly uniform since the jump of the modulus 
of the magnetisation during the reversal is about 
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Fig. 7. Time dependence in 10^ MCS of the modulus and 
the Mx, My components of the magnetisation per spin at 
ksT/J = 0.01 (left side) and fesT/J = 0.1 (right side) 
(N = 33 and D = J). 

10 % at ksT/J = 0.01 (for D = 0.2 J, the jump is 
around 6 %). The reversal duration for D = 0.3 J 
at ksT/J = 0.01 is about 2 x 10* MCS which is of 
the same order of magnitude as for £) = 0.1 J and 
0.2 J. 

Fori? = J {R/5 ~ 0.9) the material of our model 
is harder than Co (the quality factor is about 0.8). 
As it can be seen in Fig. 7, the reversal is clearly 
not uniform (the jump is about 36 %). Because 
of strong anisotropy, the spins behave roughly as 
Ising spins and the magnetisation is always along 
the X axis except short time scale fluctuations. This 
was not the case for D = 0.1 J (see Fig. 2). Ac- 
tually the spins near the surface of the nanoparti- 
cle firstly flip whereas the other spins with higher 
coordination number flip later, that is for a larger 
value of the field. At ksT/J = 0.01, the magneti- 
sation reversal occurs in only two steps indicating 
that the nanoparticle behaves as a core surrounded 
by a surface shell. At fc^T/J = 0.1, the reversal is 
more continuous due to thermal activation. 

To check the vahdity of Eq. (6) for D = 0.2 J 
and 0.3 J, we have plotted the mean switching field 
values using the same procedure as in the previ- 
ous sections. Although the reversal is not perfectly 
uniform, it is possible to put the data points cor- 
responding to ?,'/,/< 10^^ on a single straight line 
(Fig. 8). As for D = 0.1 J, the scaling behaviour 
is satisfied for all temperatures. The results of the 



Fig. 8. Scaling plot of the mean switching field val- 
ues of a 33-spin particle with D = 0.1 J for 
5 X IQ-" < v/J < 10-^, D = 0.2 J and 0.3 J for 
5 X 10-10 < v/J < 10"'^ and 0.01 < fcsT/J < 0.1. 

Table 3 

Parameters of the fits and estimates of the switching field 
at /sT and the energy barrier in zero field for a 33-spin 
particle with different values of the anisotropy constant. 



D/J 


0.1 


0.2 


0.3 




0.215 


0.430 


0.645 


c/ks 


10--^ 


10"^ 


10~2 


axis intercept 


0.219 


0.432 


0.646 


slope 


-0.132 


-0.207 


-0.252 


AEo/J 


2.8 


4.3 


6.6 


To 


40 


50 


5 



two fits are reported in TABLE 3 in comparison 
with the case D = 0.1 J. As for D = 0.1 J, 
a self-consistent fit is obtained with H^.^^ slightly 
higher than the expected value but our estimates of 
Hgyj are proportional to the anisotropy constant in 
agreement with the Stoner-Wohlfarth model. The 
extracted AEq is lower than the expected value 
again, the increase of the misfit with the anisotropy 
constant being attributed to the evolution towards 
a non-uniform rotation. 

For D = J, it is not possible to put the low tem- 
perature data points on a single straight line. The 
best fit is shown in Fig. 9. Note that the switching 
field values are significantly lower than that pre- 
dicted in the case of uniform rotation {H^^ /J = 2) . 
Actually, at very low temperature, only spin flips 
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Fig. 9. Scaling plot of the mean switching field values of a 
33-spin particle with D = J for 5 x lO"!" < vjJ < lO"*^ 
and 0.01 < fesT/J < 0.1. 

can occur, so the switching field does not depend 
anymore on the anisotropy constant but rather 
on the exchange interaction and the coordination 
number. We expect that the numerical estimate of 
the switching field at very low temperature would 
have not been changed increasing the value of the 
anisotropy constant. 



5. Conclusion 

Using a MC single spin rotation algorithm, we 
have realized the Neel-Brown model of the mag- 
netisation reversal for nanoparticles with realistic 
values of R/X and R/S. By choosing the temper- 
ature and the field sweeping rate appropriately, it 
has been possible to check the validity of the ap- 
proximated expression of the mean switching field. 
For reasonable anisotropy (£> < 0.3 J), the scaling 
behaviour is satisfied at very small field sweeping 
rates {v/J < 10~*) for all temperatures studied 
except the case N = 7 and D = 0.1 J for which 
the superparamagnetic regime can be observed. On 
the other hand, for high anisotropic nanoparticles 
{D = J), the magnetisation reversal is not uniform 
and the mean switching field values can not be fit- 
ted anymore following the proposed approximated 
expression. 

Let us briefiy remind that our estimates of H^^ 



are independent of the size and vary linearly with 
the anisotropy constant, as expected. The discrep- 
ancy between the estimate of AEq and its expected 
value in all cases might be due to the spin rotation 
procedure. 

In a near future, it is planned to investigate 
the effect of physical parameters, such as surface 
anisotropy and dipolar interactions, in relation 
with the shape of the particle. However, one has 
to be aware of the very large increase of CPU time 
by taking into account the dipolar interactions. 
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